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ON  TIE  CONCEPT  OF 
STABILITY  OF  INELASTIC  SYSTEMS 
By  D.  C.  Drucker^-  and  E«  T»  Onat^ 


Abstract 

Simple  models  are  employed  to  bring  out  the  large  and  im- 
portant differences  between  buckling  in  the  plastic  range  and  classical 
elastic  instability.  Static  and  kinetic  criteria  are  compared  and 
their  inter-relation  discussed.  Non-linear  behavior  in  particular 
is  oft^n  found  to  be  the  key  to  the  physically  valid  solution.  The 
non-conse  vative  nature  of  plastic  deformation  in  itself  or  in  com- 
bination with  the  non-linearity  requires  concepts  not  found  in  clas- 
sical approaches.  Conversely,  the  classical  linearized  condition  of 
neutral  equilibrium  Is  really  not  relevant  in  inelastic  buckling. 

Plastic  buckling  loads  arc  not  uniquely  defined  but  cover  a range  of 
values  and  are  often  more  properly  thought  of  as  maximum  loads  for 
some  reasonable  Initial  imperfection  in  geometry  or  dynamic  disturbance. 

The  models  indicate  that  basically  the  same  information  is 
obtained  from  essentially  stetic  systems  by  assuming  initial  imperfec- 
tion in  geometric  form  a;  by  assuming  dynamic  disturbances.  One 
aporoach  complements  the  other  and  Noth  are  helpful  in  obtaining  an 
understanding  of  the  physical  phenomena. 
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Introduction 

In  the  analysis  of  structures,  as  in  most  branches  of 
engineering  and  science,  problems  are  solved  by  simplifying  then 
enormously.  The  geometry  of  the  structure,  the  material  of  which 
it  is  composed,  and  the  loads  applied  are  all  strongly  idealized. 

A rather  extreme  example  is  a riveted  truss  which  is  often  analyzed 
as  pin-connected  and  as  though  all  axial  forces  in  the  members  were 
applied  along  the  assumed  straight  centerline  of  each  and  were  uni- 
formly distributed  over  the  cross-section.  Justification  of  any 
idealization  requires  an  analysis  of  the  real  system  and  an  inter- 
pretation of  the  results.  Should  reasonable  deviations  produce  large 
changes  in  the  result  the  idealization  is  not  permissible.  Clearly, 
in  the  example  just  cited  the  truss  members  will  have  appreciable 
end  moments  so  that  the  stress  will  not  be  distributed  uniformly. 

Largo  stress  concentrations  will  occur  at  rivet  holes  and  other 
discontinuities.  In  a strict  sense  the  idealizations  are  entirely 
in  error.  Fro  ’ the  practical  point  of  view,  hov/ever,  if  allowable 
stresses  arc  based  on  experience  as  codified  in  engineering  snecifica- 
tions,  it  will  often  bo  "ound  that  the  simplified  analysis  is  adequate. 

The  basic  clastic  buc'/ling  problem  of  the  Euler  strut  is 
more  relevant  to  the  present  discussion.  The  assumptions  made  are 
that  the  column  is  perfectly  straight,  that  the  load  is  static  and 
is  applied  along  the  centerline,  and  that  the  material  is  homogeneous, 
linearly  clastic,  and  free  of  initial  stress.  Computation  of  critical 
load  is  made  v/ith  a linearized  instead  of  exact  expression  for  curva- 
ture while  a condition  of  neutral  equilibrium  is  sought  under  constant 


load.  Neutral  signifies  that  equilibrium  is  possible  in  neighboring 
deformed  positions  as  well  as  in  the  ideal  configuration.  An  equiva- 
lent approach  is  to  compute  the  static  force  or  to  compute  the  impulse 
needed  to  displace  the  strut  laterally  and  then  to  define  the  critical 
load  as  the  load  at  which  zero  disturbance  is  needed. 

The  terms  instability  and  buckling  load  can  thus  be  given 
precise  meaning  for  an  ideal  Euler  elastic  strut  without  any  con- 
sideration of  the  behavior  of  real  columns.  However,  an  engineer 
of  skeptical  turn  of  mind  encountering  such  a calculation  for  the 
first  tine  could  well  be  excused  if  he  ignored  the  results  completely. 
Indeed  they  were  generally  ignored  until  aircraft  sections  emphasized 
their  significance.  There  is  no  obvious  or  intuitive  reason  for 
accepting  the  validity  of  the  idealizations.  Rolled  and  extruded 
sections  have  appreciable  initial  stress  and  are  imperfect  in  shape 
and  form.  Loads  are  not  absolutely  axial  nor  are  they  ordinarily 
constant  in  magnitude. 

Fortunately,  in  this  simple  problem,  idealizations  and 
imperfections  can  be  studied  analytically.  When  Euler  solved  the 
problem  originally,  he  used  the  exact  expression  for  curvature  so 
that  its  linearization  is  known  to  be  permissible.  Also,  the  eccen- 
trically loaded,  imperfect  column  can  be  treated.  The  usefulness  of 
the  Euler  critical  load  computation  lies  in  the  fact  that  the  calcula- 
tions based  on  real  columns  show  the  critical  load  to  be  a reasonable 
limiting  load  as  long  as  the  stresses  induced  are  below  the  elastic 
limit  and  the  imperfections  are  small,  Fig.  1. 
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Similar  types  of  calculations  have  proved  equally  successful 
for  other  structural  elements  where  small  deviations  from  the  ideal- 
ized condition  also  do  not  produce  large  differences  in  behavior.  In 
a way,  this  successful  analysis  of  linearized,  idealized,  and  simpli- 
fied systems  was  unfortunate.  It  seems  to  have  led  to  the  idea  that 
the  classical  linearized  theory  with  its  static  criterion  of  neutral 
equilibrium  solved  the  real  problem.  As  an  actual  imperfect  system 
is  so  troublesome,  even  in  the  elastic  range,  it  is  easy  to  see  how 
such  a situation  could  develop.  When  low  values  were  obtained  experi- 
mentally for  the  buckling  loads  of  cylindrical  and  spherical  shells 
the  warning  was  not  taven  seriously  in  general.  Donnell  did  attribute 
the  discrepancy  to  initial  imperfection  but  harman  and  Tsien*  explained 
the  result  with  a large  deflection  theory,  Fig.  2.  Although  the  two 
approaches  are  comparable  in  some  ways,  it  was  the  buckling  computa- 
tion which  had  the  more  popular  appeal. 

The  picture  began  to  change  when  Shanley**  introduced  the 
concept  of  considering  the  loading  process  itself  by  returning  to.  the 
strut  problem  and  following  what  happens  as  the  load  is  increased  and 
the  elastic  limit  is  exceeded.  He  demonstrated  conclusively  that  the 
classical  buckling  approach,  which  gives  the  critical  load  P^.  is  not 
appropriate  for  the  plastic  range.  Perhaps,  however,-  too  much  atten- 
tion was  paid  to  the  remarkable  proof  that  an  initially  perfect 
column  could  start  to  bend  at  the  tangent  modulus  load,  P^ . The 

*Th.  V.  Kerman  and  H.  Tsien,  Jour. Aero. Sci.  Vol  7,  PP  43-50,  1939 

**F.R.  Shanley,  Jour. Aoro. Sci.  Vol  13,  p.  678,  1946 

Jour. Aero. Sci.  Vol  l4,  p.  26l,  194-7 
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buckling  aspects  of  Shanley's  problem  seemed  to  overshadow  the  max- 
imum load  computation,  or  more  generally,  his  load-reflection  relation, _ 
possibly  because  the  spread  between  and  is  so  small  for  a column. 
As  already  mentioned,  initial  imperfections  are  always  present  so  that 
basically  the  question  of  practical  importance  is  how  the  deflection 
grows  with  load  and  what  is  the  maximum  load,  which  can  be  carried. 
Pearson*  contributed  to  the  overall  picture  by  proving  that  Shanley's 
load-deflection  curve  for  an  initially  straight  column  is  the  limiting 
load-deflection  curve  for  an  imperfect  column  as  the  deviation  from 
straightness  approaches  zero,  Fig.  3. 

Onat  and  Drucker**  added  the  more  elaborate  example  of  a 
cruciform  column  in  the  plastic  range  which  fails  by  twisting,  Fig.  4-. 
This  plate  problem  solved  by  small  deflection  theory  demonstrated  the 
need  to -take  into  account  what  might  bo  termod  small  but  finite  de- 
formations which  occur  as  the  loading  proceeds  if  extremely  small 
initial  imperfections  arc  present.  Kero  the  spread  between  P;vj  and  P^ 
is  large  and  of  real  practical  significance  similar  in  a \ ly  to  the 
shell  problem  of  Fig.  2. 

Ziegler***  focussed  attention  on  dynamic  loading  or  dynamic 
disturbance  superposed  on  static  loading  and  pointed  up  the  shortcomings 


*C.E.  Pearson,  Jour. Aero. Sci.  Vol  17,  p.  4-17,  1950 

**E.T,  Onat  and  D.C.  Drucker,  Jour. Aero. Sci.  Vol  20,  pp.l8l-l86,  1953 

***H,  Ziegler,  Ingenieur-Arc^iv,  Vol  20.  pp. 4-9-56,  1952 

Ziegler,  Zeitschrift  fur  angewandte  Mathematik  and  Physik,  Vol  4-, 
Fasc.  2 and  3,  1953 
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of  classical  theory.  A kinetic  criterion  of  instability  can  be 
employed  to  obtain  correct  answers  for  systems  with  perfect  geometry. 
The  necessity  or  convenience  of  a kinetic  criterion  for  essentially 
static  problems,  however,  remains  to  be  considered. 

To  summarize  these  introductory  remarks,  the  real  problem 
of  instability  involves! 

1)  imperfection  of  geometry  of  structure 

2)  imperfection  of  loading 

3)  dynamic  disturbances 

L:-)  inhomogeneity,  residual  stresses,  and  similar 
ever-present  imperfections 

In  general  calculations  must  be  made  taking  into  account  the  effect 
of  small  but  finite  deformations  and  possibly  large  deformation  as 
well  in  a fev  cases.  On  the  other  hand,  classical  linearized  theory 
for  the  idealized  system  computes  instability  for  infinitesimal  de- 
formation or  dynamic  distrubance  under  constant  load.  If  the  loading 
is  conservative  and  the  system  is  linearly  elastic,  such  an  analysis 
is  often  satisfactory.  The  same  critical  load  is  obtained  for  both 
infinitesimal  deformation  and  dynamic  disturbance.  As  inelastic 
systems  are  both  non-linear  and  non-conservative  (path  dependent),  it 
is  not  reasonable  to  expect  the  classical  linear  type  of  theory  to 
produce  significant  results. 

In  the  Shanley  example  Fig.  3»  it  is  non-linearity  which  is 
responsible  for  deviation  from  straightness  at  the  tangent  modulus 
load,  non-conservatism  which  produces  Pj^.  Non-linearity  is  so  promin- 
ent in  the  elastic-plastic  range  because  a relatively  small  amount  of 
deformation  produces  the  change  from  elastic  response  with  its  high 


I 


I 

t 

! 
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modulus  of  elasticity  to  the  plastic  resnonse  at  low  modulus. 

The  meaning  of  instability  is  itself  often  obscure  and  a 
matter  of  new  definition.  Loads  betv/een  P^,  and  Pj^,  Fig.  3 are  stable 
in  one  sense  unstable  in  another.  Dynamic  disturbances  and  geometric 
imperfections  must  be  considered  and  the  history  of  the  deformation 
must  be  follov/ed. 


Non-Linearity,  Static  Analysis 


As  stated  plastic  action  is  always  both  non-linear  and  non- 
conservative. The  first  i odel  to  be  considered,  rig.  has  both 
characteristics  but  it  i ; non-linearity  which,  provides  the  interesting 
features  of  its  behavior.  Rod  CA  is  rigid,  the  pin  at  0 is  friction- 
less  and  the  spring  attached  at  A lias  an  elastic  and  a elastic  range 
of  force  v shown  schematically  by  the  full  line  in  Fig.  6.  Small 
displacements  only  will  be  analyzed  so  that  1,0  replaces  L sin  0, 

L cos  0 is  taken  as  I.  for  the  lever  arm  of  F,  and  t!  e position  of  A 

O 

belov;  its  maximum  possible  height  is  approximated  by  L0  /2.  Upon 
unloading,  the  force-displacement  curve  is  straight  and  parallel  to 
the  original  elastic  line  as  shov.-n.  The  actual  force  displacement 
curve  'rill  be  replaced  by  the  broken  line  segments  I3C,  CD,  DD*  , etc. 
for  convenience  of  description  and.  of  algebraic  manipulation. 

The  mathematical  expression  of  the  simplified  force-displace- 
ment relation  depends  uoon  wnether  the  material  is  behaving  elastical- 
ly or  plastically  and  is 

dF  = 1:  LdO 


(1) 


Expression  (1)  applies  for  FdF  > 0 and  0 >0^  on  first  loading  or  more 
generally  FdF  > 0 and  I w I s Fy  and  greater  than  any  value  reached 
previously.  Expression  (2)  applied  for  FdF < C or  for  I F I < F or  any 
previously  attained  value. 

Suppose  the t the  spring  is  so  adjusted  that  F = 0 when 
0=0.  The  bar  OA  will  then  he  in  equilibrium  in  the  vertical  posi- 
tion for  all  values  of  load,  P.  The  question  of  tine  stability  of  the 
equilibrium  is  not  as  trivial  as  night  at  first  appear.  If  the  usual 

static  criterion  is  employed,  OA  is  rotated  through  a very  small 

2 

angle  0 and  the  work  done  by  F,  PL©  /2  is  compared  v.ith  the  energy 
stored  by  the  spring,  FLO/2  where  F = kLG.  Equating  the  two  terms 
leads  to  the  critical  load 

PK  = kL  (3) 


The  Shanley  concept  of  increasing  P as  0 is  introduced  does  not  change 
the  result  a^  ail.  It  is  tie  stability  of  deflected  Dositions  which 
is  the  ’ey  to  the  physically  significant  behavior  of  t.he  system. 
Equilibrium  squires  PL©  = FI.  or 

F = F©  (4) 

For  t elastic  range,  ©<©-,,  F = kLO  and  P = kL.  When  © reaches  or 
exceeds  ©j_,  F = Fy  + ktL(©  - 0^)  whore  Fy  = kL©^. 


•Hv+ur*  imm 
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F6  = kL©x  + ktL(©  - ©1)  = (k  - kt)L©1  + ktL© 


= kL©  - (k  - k+)  L (©  - ©,) 


(5) 


as  depicted  in  Fig.  7 by  lines  ©0  = 0. 

If  the  spring  is  so  adjusted  that  F = 0 at  © = ©0, 

F = kL( © - ©0)  = P© 


and 


e = rz 

1 kL 


l-£- 

:K 


(6) 


in  the  elastic  range,  Icy end  Fy,  where  as  shown  dashed  in  Fig.  7a 

P 


pk  ©o+ei 


F - kL©i  + ktL  (©  - ©D  - ©x)  r-  P« 


so  that  to  maintain  equilibrium  F would  have  to  decrease  as  © in- 
creases. The  maximum  value  of  P therefore  is  given  by 


M 


Qo+el 


(7) 


with  the  restriction  PN>  > k^L. 

In  a sense  Figs.  7(a)  and  7(b) 


rhich  differ  only  in 


, 


horizontal  scale  are  a graphical  representation  of  the  two  extreme 
possibilities.  If  ©^  is  large  compared  v/j  th  probable  initial  and 
subsequent  deflections  P = kL  is  a valid  buckling  load  and  will  be 
found  errneri-ner.kc.-lly.  On.  the  contrary,  if  ©^  is  very  small  compared 
with  probable  initial  deflection  or  equivalent  eccentricity  or 
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Inclination  of  load  P,  the  experimentally  determined  critical  load 
or  the  one  found  in  practice  will  be  close  to  k^L.  When  k^.<<k, 
the  difference  between  the  extremes  is  very  large.  To  this  extent, 
the  physical  phenomenon  may  be  said  to  be  strongly  dependent  upon 
initial  imperfection.  In  particular  cases  as  for  the  cruciform, Fig.  *+, 
it  may  turn  out  that  the  maximum  load  is  bounded  quite  closely  in 
practice. 

If  the  model  is  taken  as  important  in  itself  instead  of 
simply  a schematic  illustration  it  is  interesting  to  see  how  the 
extreme  cases  arise.  Fy  -=  kL©^  may  be  rewritten  as 


where  6e  is  the  maximum  elastic  elongation  of  the  spring.  A long  soft 
spring  will  provide  a S>e  of  several  inches  and  P,,  will  be  obtained. 
Short  snubbing  bars  or  wires,  on  the  other  hand,  permit  a t of  the 
order  of  thousandths  of  an  inch  and  unavoidable  measures  in  the  bar 
and  the  loading  will  give  an  effective  ©Q  many  times  ©^.  Pj.;  will 
then  tend  to  be  close  to  k^L. 


Although  Fig.,  6 shows  the  unloading  behavior  to  be  non- 
conservative, nothing  in  this  section  involved  unloading  so  that  a 


non-linear  elastic  material  would  exhibit  exactly  the  same  effects. 


Non-Linearity,  Fine tic  Analysis 

The  model  of  Fig.  5 will  now  be  considered  from  the  alterna- 
tive kinetic  point  of  view.  Now  the  question  to  bo  asked  is  whether 
a dynamic  disturbance  applied  to  a perfectly  aligned  system  in 
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equilibrium  produces  bounded  or  unbounded  oscillations  within  the  frame 
work  of  small  deflection  theory.  The  previous  static  analysis  makes 
it  quite  clear  that  the  answer  will  again  depend  upon  the  magnitude 
of.  the  disturbance  just  as  the  size  of  ©0  was  significant. 

Suppose  that  under  constant  load  P,  with  © = 0,  the  bar 
OA  is  given  an  initial  velocity  ©0.  Conservation  of  energy  gives, 
for  small  ©, 

I(62  - e02)  = PL©2/2  - kL©2/2  (9) 

for  © < ©2 
and 

I(©2  - ©02)  = PL©2/2  - [ kL©12/2  + kLQ1  (©-©-,) 

+ ktL  (©  - q±)2/2  ] (10) 

for  0 > ©^  and  ©_>©•  I is  the  moment  of  inertia  of  the  bar  OA  about 

0. 

The  system  is  unstable,  that  is  motion  in  the  initial 
direction  will  not  stop,  0 > 0,  if  the  energy  stored  and  dissipated 
in  the  spring  can  not  be  as  large  as 

I©02  + PL©2/2  (11) 

Obviously,  from  (Q)  and  (10)  ? < Pj.'  = kL  Ls  a requirement  for  stability 
At  P = p£,  the  slightest  ©D  will  cause  collapse,  at  P > P^  the  system 
runs  away  quickly.  A plot  of  © is  © is  most  convenient  for  exhibiting 
such  properties,  Fig.  8. 

• p 

If  P <kj.L,  there  is  no  0O  which  can  not  be  absorbed  in 
accordance  with  small  deflection  Equation  (10).  The  system  is 
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therefore  stable  for  F<ktL  and  the  0,  6 diagram  is  composed  of 
portions  of  ellipses,  Fig.  9.  For  ktL  < P < kL  there  will  be  levels 
of  initial  disturbance  above  which  the  system  will  not  recover.  Below 
those  values  of  Q0  the  bar  will  come  to  rest  and  return.  It  is  clear 
on  physical  grounds  that  on  the  return  motion  the  system  will  again 
stop  and  oscillate  back  and  forth.  Unless  the  system  collapses  on 
the  first  try  it  is  stable,  Fig.  10.  It  is  for  this  reason  that  on 
this  and  analogous  cases  whon  tho  external  forces  applied  are  con- 
servative, static  and  kinetic  analyses  are  equivalent. 

Fig.  10  also  s'-ows  the  effect  of  combinations  of  nitial 
velocity  and  initial  dis ^lacenent . 

NTon-Conservativc  Aspects 

Fig.  11  shows  the  model  employed  by  Duberg  and  Wilder*1,  fol- 
lowing closely  the  model  and  concept  of  Stanley,  to  investigate  the 
effect  of  initial  imperfection  in  columns.  DLGR  is  rigid  and  G is 
constrained  to  more  vertically  only  by  means  of  a frictionloss  guide 
so  that  the  system  h--s  two  degrees  of  freedom  y,  9 . L and  B rest 
on  two  identical  short  brrs  of  work-hardening  material  whose  force- 
displacement  diagrams  are  each,  idealized  as  two  straight  lines,  Fig. 
12.  If  a bar  bias  been  loaded  to  B and  loading  continues 

F = Fr  + kt  ( 6 - 63)  (12) 

while  if  load  is  remov  ■••1 


*J.E.  Duberg  and  "ilier,  Jour. A ro.Gci.  Vol  17,  ?,?23,  195' 


utfLa*  SHtuik 


■vvh 
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F=FQ-k(6B-6) 


The  static  analysis  is  well  known  and  basically  has  already 


been  described  in  Fig.  3 where  6 should  be  replaced  by  ©„  The  load 


rotation  relation  for  the  simplified  model  is  given  in  Fig.  13.  As 


initially  imperfect  system,  © = © , will  rotate  more  as  load  is 


applied  and  will  have  deviated  appreciably  fron  the  vertical  by  the 


time  the  tangent  modulus  load,  P^,,  Fig.  14,  is  reached, 


pt  = ££ 


(14) 


A nerfect  syste*  can  deviate  from  the  upright  position  at  P-j.  and  the 
slope  of  the  P,  © line  will  be  progressively  flatter  the  larger  the 


value  of  P at  which  rotation  is  permitted  (see  line  segments  between 


PT  and  PK  of  Fig.  3).  At  the  reduced  modulus  load  P = PK,  Fig.  1*+, 


the  slope  would  be  zero 


k*b“  2k 


pk  = -l~  irr+K 


(15) 


The  third  critical  load  of  interest  is  the  Fulcr  clastic 


buckling  lo^d  which  can  be  obtained  by  putting  = k in  (1*+)  or  (15) 


p _ kb£ 


(16) 


The  rr-  ceding  discussion  has  implicitly  assumed  that  F-j.  and 
therefore  P^  and  are  above  the  yield  force  PF^,  Fig.  12.  Should 
P^,  be  loss  than  2F  , the  initial  imperfections  may  be  small  enough 


so  that  the  material  may  not  bo  made  aware  of  the  plastic  deformation 


it  could  experience.  Nothing  unusual  will  then  happen  at  the  tangent 
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modulus  load.  The  obvious  example  is  a perfectly  plastic  or  non-work- 
hardening material,  kt  = 0,  for  which  P-p  = 0.  In  what  follows,  t5  ore- 
fere,  P<|>  > ~V  y* 

A dynamic  investigation  of  the  geometrically  perfect  model 
I3  more  interesting  than  the  well-explored  static  study.  As  in  the 
analysis  of  the  first  model,  the  consideration  of  a kinetic  criterion 
of  stability  loads  to  investigation  of  the  motion  which  follows 
arbitrary  initial  disturbances,  and  the  question  to  be  answered  is 
whether  or  not  the  ensuing  motion  is  limited.  The  system  has  two 
degrees  of  fr  iodon,  one  of  motion  of  0 vertically,  y,  and  one  of 

rotation,  ©,  and  the  arbitrary  initial  dynamic  disturbance  will  be 

• • 

characterized  by  yQ,  ©c  where  dot  indicates  time  derivative.  TMs 
pair  of  values  determines  the  initial  instantaneous  center  of  rotation 
C and  of  course  the  initial  angular  velocity  cf  the  rigid  bar  DLhR. 

If  the  x-coordir.ate  cf  C is  between  L and  R,  one  supporting  bar  will 
start  to  unload  and  the  other  to  load.  On  the  oth^r  hand,  if  the 
x-coordinat a of  C lies  outside  L to  R both  bars  will  start  to  lead  or 
unload  depending  upon  the  sign  of  ©_.  I 

Designating  the  kinetic  energy  of  the  system  by  T,  for 
small  displacements  the  change  in  kinetic  energy  is 

2 

AT  = T - T0  = P(y  + LG-)  + ?L  VL  dt  + 

where  F and  V are  the  instantaneous  values  of  force  and  velocity  at 
L and  R.  Denoting  the  bar  constants  by  k^  and  kp,  each  of  which  may 
be  either  k or  k^.,  and  considering  the  very  beginning  of  the  mot  Lon, 
first  order  terms  cancel  out  and 


T = ^.|(y^e)  - ^ (y  - b 0)2# 


(18) 


The  model  therefore  behaves  like  a linear  elastic  system  at  the  begin- 
ning of  the  motion  except  that  the  non-conservative  stress-strain 
relations  determine  the  proper  value  of  tho  k's.  A decrease  in  kinetic 
energy,  AT  < 0,  for  all  possible  choices  of  y and  0 indicates  that  the 
system  is  probably  stable.  An  increase  does  not  necessarily  indicate 
true  instability  because  even  within  the  limitation  of  small  deflection 
theory  non-conservatistt  means  that  (l8)  need  not  hold  as  the  motion 
proceeds . 


If  both  support  bars  load,  k^  = ^ and  y2^  There- 

fore from  (18),  At  will  bo  negative  (indicating  stability)  for 
k«-b2 

P < PT  = — . For  P > PT  there  will  be  a set  of  initial  disturbances 

for  which  the  kinetic  energy  of  the  system  will  increase  at  the  begin- 
ning of  the  motion. 

If  0Q  is  positive  and  one  support  bar  unloads  while  ore 

other  increases  its  load,  k^  - k^.,  k^  = k,  and  y < ^ ©.  With  these 

conditions  Equation  (1°)  shows  that  AT  will  be  negative  for 
ktb2  2k 

P < Fv  = — £ — r r • -:’he  1 inetic  energy  of  a disturbance  of  this  type 

K L kt+k 

will  not  grow  until  the  load  exceeds  tho  reduced  modulus  load.  What 
is  of  possibly  greater  significance  here  is  that  for  P>PIc,At>0  and 
the  model  will  collapse  because  the  k' s do  not  change  as  the  incision 
proceeds. 

If  the  artifical  case  is  chosen  in  which  both  bars  unload 
= k and  y < - ^ . A?  will  be  negative  for 
slowing  down  does  not  necessarily  indicate 


initially,  k^ 


P <Px 


kb£ 


■ 

Sue 


All-92 


16 


stability  in  this  case.  Obviously,  after  reversing  its  motion,  the 
system  will  collapse  if  P >PK. 

In  general,  therefore  the  complete  motion  of  the  system 
must  be  examined  under  all  possible  disturbances  to  understand  the 
significance  of  the  critical  loads  obtained.  As  illustrated  in  Fig. 15 
because  k^.  is  constant  the  motion  sottles  down  when  P-p  £ P_£P^ 
whether  AT  is  initially  positive  or  negative.  With  the  usual  curved 
stress-strain  diagram  it  is  clear  that  collapse  would  occur  for  P 

loss  than  P^  by  an  amount  depending  upon  the  magnitude  of  the  initial 

• • 

values  ©0  , yG.  Nevcrtaeiess  it  is  interesting  and  important  to  note 
that  all  three  critical  loads,  P^,  P^,  and  Fp  can  be  obtained  from 
an  analysis  of  a geometrically  perfect  system  subjected  to  infinites- 
imal dynamic  disturbance. 
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